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ABSTRACT 

Cosmological observations from ground at millimetre and sub-millimetre wavelengths 
are affected by atmospheric absorption and consequent emission. The low and high 
frequency (sky noise) fluctuations of atmospheric performance imply careful observa- 
tional strategies and/or instrument technical solutions. Measurements of atmospheric 
emission spectra are necessary for accurate calibration procedures as well as for site 
testing statistics. CASPER2, an instrument to explore the 90 -j- 450 GHz (3-^15 cm -1 ) 
spectral region, was developed and verified its operation in the Alps. A Martin-Puplett 
Interferometer (MPI) operates comparing sky radiation, coming from a field of view 
(fov) of 28 arcminutes (FWHM) collected by a 62-cm in diameter Pressman-Camichel 
telescope, with a reference source. The two output ports of the interferometer are de- 
tected by two bolometers cooled down to 300 mK inside a wet cryostat. Three differ- 
ent and complementary interferometric techniques can be performed with CASPER2: 
Amplitude Modulation (AM), Fast-Scan (FS) and Phase Modulation (PM). An al- 
tazimuthal mount allows the sky pointing, possibly co-alligned to the optical axis of 
the 2.6-m in diameter telescope of MITO (Millimetre and Infrared Testagrigia Ob- 
servatory, Italy). Optimal timescale to average acquired spectra is inferred by Allan 
variance analysis at 5 fiducial frequencies. We present the motivation for and design 
of the atmospheric spectrometer CASPER2. The adopted procedure to calibrate the 
instrument and preliminary performance of the instrument are described. Instrument 
capabilities were checked during the summer observational campaign at MITO in July 
2010 by measuring atmospheric emission spectra with the three different procedures. 

Key words: Instrumentation: interferometers - Cosmology: observations - Site test- 
ing 



Ground-based cosmological observations can be carried 
out in the millimetre and sub-millimetre bands (hereafter 
mm and sub-mm) but need a continuous monitoring of the 
atmospheric contribution. Accuracy in the calibration pro- 
cedure towards known photometric sources results from a 
good knowledge of transmission and emission along this op- 
tical path. Sky brightness can be suitably monitored in this 
frequency range by spectrometers designed only with this 
goal (see as exampl e the FTS interferometers employed in 
iMatsuo et all l|l998h . iMatsushita et al.l l| 19991 ') . iPaine et all 
(2000)) or, more generally, to spectral observations of sky 
objects when optically matched to existing telescopes (see 
the desc ription of the instrument inst alled at the CSO focal 
plane in ISerabvn fc Weissteinl (|l996l )h We have developed 
an instrument, CASPER2, to record atmospheric emission 



spectra in the millimetre band for assisting cosmological ob- 
servations with the 2.6-m in diameter Cass egrain telescope 
at MI TO (Testa Grigia, Italy, 3480 m a.s.l.) |De Petris et all 
( |2007ft ). It is conceptually similar to CASPER, the experi- 
me nt proposed for t he ita lian-french Antarctic base, Dome 
C (|De Petris etafl (120051)), but with a more limited spec- 
tral band (jDecina et al.l (|2010l 0. Such an instrument makes 
it possible to avoid specific telescope procedures, with con- 
sequent loss of observational time like skydips, and to infer 
atmosphere opacity in a wide spectral range. 
After a brief description of the importance of an instrument 
like CASPER2 in Sect. 2, the instrument concept and its 
major characteristics are described in Sect. 3. In Sect. 4 a 
discussion about the Martin-Puplett interferometer (MPI) 
and the options adopted for signal sampling are presented. 
Sect. 5 reports preliminary atmospheric spectra as recorded 
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at MITO in July 2010. Final remarks are discussed in the 
Conclusions. 



1 COSMOLOGICAL OBSERVATIONS 
REQUIREMENTS 

Ground based observations, focused on cosmological targets 
in the mm-submm bands (sub-THz frequencies), are greatly 
affected by atmospheric presence. The performance in 
this spectral region is mainly dominated by water vapour 
content modeled by a continuum-like term and by a series 
of absorption lines peaked at 183, 325, 380, 448, 557 GHz. 
Oxygen is also present as absorber at 118, 368 and 425 
GHz but, differently from H2O fluctuations, it is well-mixed 
in the atmosphere and so it mainly contributes only as 
a constant radiation background, still dependent on the 
elevation, with consequent photon noise on the detectors. 
Atmospheric synthetic spectra, derived from models, are 
a valuable support for predicting atmosphere brightness 
starting from thermodyn amic param e ter v alues, see for 
example the ATM code |Pardo et all (|200lh V However a 
validation of the model for a specific site is mandatory 
and this is possible only when on-site recorded spectra 
are available. This is what is going on at MITO with 
CASPER2. 

Short term variation of the atmospheric emission on 
different spatial scales contributes as sky noise while 
slow fluctuations of the attenuation can affect the quality 
of sky calibrations. Regarding only observations in the 
millimetre/sub-millimetre bands, two kinds of approaches 
for atmosphere monitoring are widely applied. The atmo- 
sphere is usually dealt with by an isothermal slab model 
where opacity follows a zenithal secant-law. Under this 
assumption large telescopes are periodically employed to 
perform altitude scans in the sky, i.e. skydips, to derive 
the zenith optical depth, integrated in the operating photo- 
metric band. Due to the fact that this procedure subtracts 
time from the astronomical observations, the opacity is 
alternatively derived by employing ancillary instruments 
such as tippers, GPS or water vapou r radiometers ( see for 
instance iRadford fc Holdawavl l| 19981 ). ICoster et all (|l996T ) 
and IWiedner et al.l (|200ll )). All these usually feature wide 
fields of view and limited spectral capabilities, unmatched 
to those of scientific instruments. This makes the assessment 
of the effective site quality (i.e. in-band transparency and 
fluctuations at the angular scales of interest) a non-trivial 
task. A spectrometer is the natural solution for a continuous 
sky monitoring over a wide spectral band and an MPI is an 
efficient solution for reaching this goal. Furthermore if the 
spectrometer is optically matched to an its own medium size 
telescope it is possible to observe the sky towards a specific 
direction and at an angular scale closer to that of the 
main telescope. In any case a model for extrapolating the 
atmosphere opacity at different frequencies is mandatory. 

A semi-empiri cal procedure was proposed by 

iDe Greeori et aO (|2012) to perform an analysis of at- 



mospheric transmission employing radiosoundings data to 
feed the ATM code for generating synthetic spectra in the 
wide spectral range from 100 GHz to 2 THz. 



2 INSTRUMENT 

CASPER2 is composed of a two-mirror telescope, an inter- 
ferometer, a wet cryostat and an altazimuthal mount. The 
atmospheric emission spectra are acquired using an MPI 
that can be easily converted into a spectropolarimeter with 
many consequent scientific applications. Here the MPI is 
fully described, while in Sect. 3 the basics of the interfer- 
ometer and the calibration procedures are detailed. The de- 
tectors, two Ge-bolometers with a NEP ~ 10~ 15 W Hz 1/2 , 
are cooled down to 290 mK by a wet cryostat, using liq- 
uid nitrogen and helium combined with a He 3 fridge. The 
cryostat (I nfrared Labs, HDL-8) is identical to the one de- 
scribed in ICatalano etaD [|2004h . In Figure [t] a schematic 
CAD drawing of CASPER2 is shown while Figure [2] illus- 
trates the instrument in operation at the Testa Grigia sta- 
tion (Alps, 3480 m a.s.L). 

The optical design satisfies the expected measurement re- 
quirements: a low resolution spectrometer in the range of 
90 - 450 GHz (R ~ 50) and a small field of view (less than 
1 degree). To achieve these goals we matched a two-mirror 
telescope with an MPI. The optical layout with the main 
components, as labelled in the text, is shown in Figure [3] 
A f/3.5 Pressman-Camichel telescope collects sky radiation 
and feeds one of the two input ports of the MPI. With this 
optical solution the wavefront aberrations introduced by the 
62-cm concave spherical primary mirror (Ml) are compen- 
sated by a 12-cm convex ellipsoidal subreflector (M2). The 
primary mirror is underilluminated by the secondary mirror, 
which operates as an aperture stop of the telescope alone 
resulting in an entrance pupil 46-cm in diameter. Both the 
mirrors were manufactured in an aluminum alloy, ensuring 
a fast and homogeneous thermalization of the surfaces. The 
mirrors were carefully polished to reflect visible light: the 
final rms surface roughness is lower than 0.1 fim. The pri- 
mary mirror has been manufactured by Officine Ottico Mec- 
caniche Marcon di San Dona di Piave (Italy) while the el- 
lipsoidal subreflector by the INFN machine workshop at the 
Department of Physics in Rome. The secondary mirror is 
maintained in the right position along the telescope axis by 
a 2 cm thick polystyrene foam plate (ST) (BASF, Styrodur 
3035N), allowing a null obscuration due to the subreflector 
support also avoiding the possible consequent diffracted ra- 
diation. Laboratory measurements show high transmission 
value (> 97 per cent) and low polarization (< 1 per cent) 
of this material in the whole band. The telescope is shielded 
from off-axis unwanted radiation (i.e. Sun during daytime 
measurements) by 8 panels with inner reflective surfaces, 0.5 
mm thick Pe raluman sheets , shap ed as vanes to operate like 
roof mirrors (|Gervasi et al.1 i|l998l )). 



2.1 Optics 

A 45 degrees tilted wire-grid (WG1) linearly polarizes the 
transmitted sky radiation focused by the telescope (sky) 
(see Fig. [3}- The radiation emitted by a reference source at 
ambient temperature (ref), realized by a disc of Eccosorb 
AN72, enters, as the second input port of the MPI, reflected 
by WG1 with a perpendicular polarization. An alternative 
colder reference source, closer to atmospheric emission, is 
under consideration. The two inputs are combined in planar 
waves, at the entrance of the MPI, by an HDPE (High Den- 
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Figure 1. CAD drawing of CASPER2. 



Table 1. CASPER2 instrumental main features. 



Telescope 


1/3.5 Pressman-Camichel 


Telescope effective focal lenght 


1621 mm 


Primary mirror diameter 


620 mm 


Primary mirror conic constant 





Primary mirror curvature radius 


978.4 mm 


Secondary mirror diameter 


120 m 


Secondary mirror conic constant 


8.86 


Secondary mirror curvature radius 


354.9 mm 


Entrance pupil diameter 


460 mm 


Field of view (FWHM) 


26 arcmin 


An 


0.05 cm 2 sr 


Interferometer 


Martin-Puplett 


Spectral range 


Channel 1: 90 -r 360 GHz 




Channel 2: 90 -f 450 GHz 


Mechanical path difference 


30 mm 


Spectral resolution 


5 GHz 


Detectors 


2 composite NTD 




bolometers @ 0.3 K 


Calibrator 


Eccosorb AN72 


Mount 


altazimuthal 


Star tracker field of view 


1 arcmin 


CCD field of view 


(14.4 x 13.6) arcmin 



the two split optical beams changing the Optical Path Dif- 
ference (OPD). RM1 is mounted on a linear stage (AICOM 
S.P.A., Model SMP-123) to perform a ±15 mm mechanical 
path difference, corresponding to a 5 GHz spectral resolu- 
tion, while RM2 is sinusoidally wobbled by a shaker (Lynge 
Dynamic System, Model 409) on a linear stage with an am- 
plitude of approximately 1 mm (see Sect 13.5.3]) . 

The radiation undergoes a polarization rotation of 90 
degrees when reflected back onto the roof mirrors. A second 
HDPE plano-convex lens (L2) focuses on the beam, exit- 
ing from the MPI, in front of the cryostat window. The last 
HDPE plano-convex lens (1/3) is mounted inside the cryo- 
stat, on the radiation shield of the helium liquid tank, cooled 
down to 1.6 K. This lens refocuses the radiation towards two 
detectors, i.e. the two output ports, after the splitting of the 
two polarization states by the third wire-grid (WG3), cooled 
down to 300 mK. We name Channel 1 the port correspond- 
ing to the transmitted radiation, and Channel 2 the other. 
All the wire-grid polarizers have 10-/im in diameter tung- 
sten wires spaced by 25-fim and show an efficiency in the 
reflected and transmitted polarizations better than 10 -3 in 
our spectral range. 

The first two lenses, LI and L2, were carefully shaped 
to image the subreflector on L3, at least for the Zero Path 
Difference (ZPD) RM1 position, defining it as the cold aper- 
ture stop of the full optical system, i.e. the Lyot stop. We 
have the possibility of selecting the last optical element as 
aperture stop only because a single pixel is present at the 
focal plane. Incidentally L3 is also the exit pupil. The opti- 
cal design was developed with ZEMAX0 

ensuring diffraction 

limited performance in the whole 90 -f- 450 GHz spectral 
range. The two bolometers are illuminated by f/3.5 Win- 
ston cones with a 10.5 mm in diameter aperture entailing a 
throughput of 0.05 cm 2 sr. 

The optical axis, from the MPI exit to the cryostat en- 
trance window, overlaps the elevation axis, ensuring a ver- 
tical position for the cryostat for all the telescope pointing 
positions. The optical consequences of this choice, as derived 
from the variable orientation between the wire grids with the 
elevation angle, are discussed in Sect. 3. 

After fixing the telescope configuration and assuming a 
fov equal to 28 arcminutes (FWHM) we limited the lower 
value for the focal length of the lens LI, /li, to satisfy the 
Jacquinot condition. The minimum spectral resolution can 
be related to telescope-MPI optical matching in the follow- 
ing way: 



5v > 



d tfs 
hi 



8 



fovft, 
hi 



(1) 



where v ma ,x is the maximum frequency (450 GHz), dtfs 
the telescope field stop diameter (13.2 mm) and f te i e the 
telescope effective focal length (1621 mm). In our case /li 
was chosen equal to 145 mm satisfying the Jacquinot condi- 
tion Up tO R~1000 at Umax- 



sity PolyEthylene) plano-convex lens (LI). A second wire- 
grid (WG2), still tilted at 45 degrees inside the MPI, is 
rotated around the optical axis of 35.26 degrees to correctly 
split the two polarization axes. It acts as beam splitter sep- 
arating the linearly polarized incoming radiations into two 
orthogonal components. Two 90 degree roof-mirrors (RM1 
and RM2), diamond-turned stainless steel, can move along 



2.2 Filters chain 

Blocking filters are employed in our system in order to re- 
duce radiative input during the different cryogenic stages 
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Figure 2. CASPER2 in operation at Testa Grigia station (3480 m 
a.s.l.) close to MITO telescope during the July 2010 observational 
campaign. 



M1 



-- sky 




RM2 



Figure 3. CASPER2 optical layout: all the labeled components 
are described in the text. 

in the cryostat (and consequently on the bolometers) while 
mesh filters to select the frequency bandwidth of interest. 
In Table [2] all the components in the filters chain are listed. 
We chose to perform measurements inside two similar band- 
widths: 90^360 GHz for Channel 1 and 90^450 GHz for 
Channel 2. Channel 2 was spectrally enlarged to explore high 
frequency atmosphere emission more prone, and so more 
sensitive, to pwv fluctuations at the expense of a larger back- 
ground emission with consequent photon noise: 0.5 nW and 
1 nW respectively, under the assumption of a pwv = 1 mm. 

The vacuum window of the cryostat is a 4 mm thick 
quartz window. The 77 K quartz filter is 3.3 mm thick with 



Table 2. Filters chain characteristics 



type 


cut (GHz) 


Temperature (K) 


Quartz window 


off (> 3000) 


300 


ARC quartz + 


off (> 1200) 


77 


black polyethylene 






Yoshinaga 


off (> 1650) 


1.6 


Mesh 


off (> 450) 


1.6 


Yoshinaga 


off (> 1500) 


0.3 


Mesh channel 1 


off (> 360) 


0.3 


Mesh channel 2 


off (> 450) 


0.3 


Winston cones 


on (< 90) 


0.3 



a diamond powder ARC. A black polyethylene filter, 0.1 mm 
thick, is employed to reduce visible and NIR background. 
At 1.6 K a Yoshinaga filter blocks fR radiation while an 
interference mesh filter limits the high frequencies at 450 
GHz. On the entrance apertures of the two Winston cones at 
300 mK, two Yoshinaga filters work as last IR blockers. Two 
final mesh filters perform the effective band selection. All the 
filters were supplied by QMC Instruments, Ltd. The lower 
limit band in frequency is dictated by the output aperture 
of the Winston cones operating as high-pass filters. 

2.3 Electronics and data acquisition 

The first stage for the read out of the bolometer signals 
are JFETs (model IFN146) in a common-drain configura- 
tion. Those unitary gain amplifiers are mounted on the He 4 
tank inside an aluminum shielded box, close to the detectors, 
but at a temperature of 120 K by self-heating. The mea- 
sured noise is ~ 3 nV/y/Hz. The expected total incident 
power on the detectors is 0.5 nW and 1 nW respectively, 
implying a sensor thermal conductivity of the order of 10 9 
W/K (Haller-Beeman). The low impedance signals from the 
JFETs feed an ambient temperature differential preampli- 
fier with a gain equal to 250 for both the channels and an 
output noise ~ 10 nV/y/Hz for frequency higher than 20 
Hz. Rechargeable batteries inside the preamplifier box sup- 
ply the bias voltage, as well as all the read out electronics 
boards and cryogenics maintenance. All the electrical con- 
nections are ensured by twisted pairs of NbTi wires (0.1 mm 
diameter) shielded within CuNi (0.03 mm thick). Depending 
on the signal modulation (see Sect. 13.5) the data acquisition 
is carried out by a low frequency (<1 Hz) data sampling 
of two synchronous demodulation amplifiers (Stanford Re- 
search fnstruments, SR850) or a high frequency (5 kHz) data 
sampling by an ADC (National Instruments PCI-6031). 

2.4 Pointing System 

The pointing system was developed with a twofold aim. The 
first is to track in the sky the field of view of the MITO tele- 
scope during its observations, exploring in this way an iden- 
tical optical path through the atmosphere, while the second 
goal is to be able to point specific directions to perform sky- 
dips for every azimuth angle. The telescope control is real- 
ized with Magellano ST7, supplied by ATEC Robotics (Ad- 
vanced Technologies for Research and Industry) . The move- 
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ment of the two axes is performed by hybrid stepper mo- 
tors (MAE, HS200) while the position is transduced by two 
incremental encoders (Baumer Electric, BHK 16OSA400). 
Both encoders ensure an accuracy of 400 steps/turn and 
they have a zero position to record and to reset the sky co- 
ordinates every observational run. This reference position 
also corresponds to the telescope rest position. 

The gear ratio for the altitude (and azimuth) axis is 
equal to 5.1 employing two gear wheels with 112 and 22 teeth 
respectively. In order to have stable and accurate move- 
ments, this low value was increased by a right angle gear 
(TLS, Sf40/PB3) with a 100:1 gear ratio. Finally the total 
gear ratio is 510. The transmission is ensured by two gear 
belts (Trasmecam, HTD 1040-8M-20 for altitude axis and 
1200-8M-20 for azimuth axis). 

A CCD camera (SBIG, ST-402ME) is used to check the 
pointing system, as star tracker. It is provided with a 135 
mm focal length lens to have a total field of view of 14.4' x 
13.6', with a field of view per pixel of about 14". The limit 
magnitude is 15. 

The co-alignment between visible and millimetre axes 
was checked in the laboratory with a fixed artificial source, 
considering its finite distance, and than verified in situ by 
pointing planets. The final error in pointing is lower than 
1 arcminute, more than enough for atmospheric emission 
measurements. 

The whole instrument is accommodated inside a 
1.6x2.9x1.4 m (WxDxH) deployable dome in PVC 
Precontract ® 502 (AMA, series 8000). 



3 THE MARTIN-PUPLETT 
INTERFEROMETER 

The MPI employs, as polarizers, three wire-grids distributed 
inside the instrument as described in Sect. 12.11 First of all 
we recall the basics of an MPI. 



3.1 MPI basics 

Stokes and Mueller matrix formalism can be used to de- 
scribe any state of polarization and superposition of beams 
when no phase relation has to be taken into account. Alter- 
natively, as for example in the analysis of an interferometer, 
one needs to perform a Jones matrix calcul ation, reverting 
to the Mueller matrix formalism afterwards (|Catalano et alj 
<|2004h l. 

We assume unpolarized radiation at the entrance of the 
two input ports and we describe it with the Stokes formalism 



Sin — Bi' t 



( 1 




V o 



(2) 



where in stands for in\ and ini, the two input sources. 
The two signals pass through WGl and are linearly 
polarized in the following way: 



S' in = M\VG\ ■ Sin = Bin 



( 1 

cos 2i?i 

sin 2$j, 
V 



(3) 



where $i n is the projected angle of the WGl principal 
axis on the plane orthogonal to the optical axis. The angles 
for the two inputs are linked as: $i„ 2 = -dim + it/2. = 
degrees corresponds to the vertical position. The Muller 
matrix used for a wire grid is 



Mwa = 



( 1 COS 2l} in 

COS 2&in COS 2 2$in 

sin 2$in cos 2i?i„ sin 2-&i n 
V 



sin 2-&in 
cos 2-&in sin 2§in 
sin 2 2#i n 




(4) 



After the MPI, each exiting beam can be described as: 



Si„ = MmPI ■ S'in = Bin 



( 1 

COS 2"&in COS 8 

- sin 2$in 
\ cos 2-din sin 8 



(5) 



where we used the Muller matrix for an ideal MPI 



A/, 



/ 1 

cos 8 
0-1 

\ sin 8 





sin 8 


- cos 8 / 



(6) 



The mechanical path difference between the two split 
beams is Aa; mec . The OPD, Ax opt , equal to 2Ax mec , is re- 
lated to the phase shift, 8, for each wavelength, A, as well 
known: 8 — 2irAx pt/X. 

WG3 splits each beam in the output ports in the fol- 
lowing way: 



So =^ 
MwG3 ■ S in = 

/ 1 + cos 2i?i„ cos 2ip cos 8 — sin 2i?m sin 2tp 

cos 2ip + cos 2"&in cos 2 2ip cos 8 — sin 2i?i„ sin 2ip cos 2ip 
sin 2ip + cos 2-&in sin 2ip cos 2ip cos 8 — sin 2#, n sin 2 2ip 





B,, 



(7) 



where o stands for out\ or out2, corresponding to the 
two output ports optically matched to Channel 1 and Chan- 
nel 2, respectively. 

Due to the fact that our detectors are only sensitive to 
signal intensity, i.e. the first element of Stokes vector, we 
can express the signals of the two output ports as: 

Iout 1 — A + +A~ cos2i?i ni cos 2<p outl cos 8— A" sin2i9 ini sin2<y3 ouil (8) 

Iout 2 = A + — A~ cos2i? ini cos 2<p outl cos5+A _ sm2-d ini sin 2<p outl (9) 

where A + = B ini + B in2 and A" = B ini — B in2 and 
rewriting ip ou t 2 = fouti + 7r/2. 

Under the assumption that $i ni = tpoutx, i.e. the polar- 
ization axis of WGl is aligned with the WG3 one, we can 
rewrite Eqs. [8] and [9] in the following way: 



hut! — A + A cos 8 
Iout 2 = A + — A~ cos 8 



(10) 
(11) 



The spectra of the two input sources are linked to Eqs. 
[8] and [9] in Sect. [531 



3.2 MPI efficiency versus telescope altitude 

The cold WG3, inside the cryostat, is at rest during tele- 
scope altitude movements. Its wires change orientation by 
an angle ipoutx with respect to the input polarizer WGl, 
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Figure 4. Dependence between ZPD output signals and altitude 
angle. The two vertical dashed lines limit the "blind" observa- 
tional range where the signal to noise ratio is lower than 3. In the 
bottom panel the normalized signals with the instrument function 
removed are shown (black circles for Channel 1 and red diamonds 
for Channel 2). 

linked to the altitude angle a as (pout = dim + a ~ f • For 
this reason we refer to the function /(a), as an instrumental 
efficiency, the signal dependence with the altitude. This has 
to be carefully taken into account to characterize the point- 
ing performance. The two outputs at ZPD, see Eqs. [5] and 
[9] can be expressed as: 

7 outl , 2 =A + ±A-/( a ) (12) 
where /(a) = cos2a. 

When a — 45 degrees the signals, even modulated by 
cos 5 along the interferogram, are null. To check this ex- 
pected instrument efficiency we filled the sky input with 
a blackbody source at ambient temperature (an Eccosorb 
AN72 sheet in front of the telescope) and, in order to get a 
high signal-to-noise ratio, we employed a Hg-lamp as refer- 
ence source instead of the disc of Eccosorb AN72. Spectra 
were acquired at altitude angles ranging from -10 to 100 de- 
grees. We also explored altitude angles lower than degrees 
and higher than 90 degrees to check the consistency and 
the repeatability of the performance. Interferograms were 
recorded in fast scan mode (see Sect. 13. 5.2[) with a 5 kHz 
scan rate, performing the time average over an acquisition 
time of about 7 minutes for altitude values ranging from 30 
degrees to 60 degrees (lower instrumental efficiency), and of 
5 minutes for the others (higher instrumental efficiency). In 
Figure [3] the two ZPD output signals are plotted versus the 
altitude angle, with error bars given by the standard devia- 
tion. At zenith, Channel 1 is proportional to A + + A - , after 
an intensity decreases down to a null signal (A - = 0), the 
rotation of polarization axis inverts the signal at horizontal 
position. The red points represent the polarized signal de- 
tected by Channel 2, as expected the trend is symmetrical to 
the other. Data points are well-fitted by the ±cos2a func- 
tion predicted by Eq. [12] (see the residuals in the bottom 
panel of FigQ] ) 

The loss of efficiency, due to the change of WGs ori- 
entations along different zenithal angles, entails an altitude 
range where the signals are negligible. 

In order to determine the width of this altitude range 
where the instrument has low efficiency, what we call the 



"blind" observational range, we estimated cv m i„ and a ma x 
corresponding to a signal-to-noise ratio ^ 3. An emission 
atmospheric spectra, as derived by ATM with a pwv = 1 
mm, was assumed as source. The vertical dot-dashed lines in 
Figure [4] limit this "blind" orientation directions: 43 degrees 
< a < 47 degrees. 

This is a restriction on the CASPER2 performance but 
it has few consequences for our purposes because the mini- 
mum altitude explored by the MfTO telescope is 42 degrees. 
This operational limit for MITO telescope is dictated by 
low atmospheric contamination requirements and by dome 
constraints. However CASPER2 can point lower altitudes, 
where the signal-to-noise ratio increases, allowing a complete 
angular range skydips. 

The expected variation in the atmospheric emission dur- 
ing skydips, due to the decreasing opacity with the altitude, 
has to be added to the instrumental efficiency. Each detector 
records a signal that can be expressed as: 

I outli2 = A+(a) ± A~a(a)/(a) (13) 

where o(a) is the atmospheric dependence and A~ is 
the input signals difference at the zenith position. We can 
notice that, due to the different spectral bands, A + and A~ 
have to be distinguished between the two channels, fn a sim- 
ple optically-thin atmospheric layer model we can assume 

a[a) = exp(— csc(cv)) (14) 

normalized at the zenith position. 

3.3 The sub-interferometer: Mickey 

A Michelson sub-interferometer, named Mickey, ensures the 
monitoring of the movements (position and velocity) of 
RM1 during fast scan measurements. The movable mir- 
ror is mounted on the back side of RM1. The source is a 
laser (Imatronic, Sigma 650/3) centered at A=650 nm. In- 
terference maxima are detected by a photodiode (Osram, 
SFH203) and then processed by a peak-counter circuit. Ac- 
curacy on the mirror position depends on the distance be- 
tween the peaks generated by Mickey, equal to A=325 nm. 

The Fourier Transform of the signal detected by Mickey 
along a return double-sided interferogram shows two max- 
ima, corresponding to a backward and forward scan velocity 
difference of about 70 /im/s. The peak-counter circuit is em- 
ployed to produce a trigger signal for data collection, fn this 
way all fast scan interferogram points are acquired at the 
same distance between each other, independently of veloc- 
ity changes due to the step motor. 

3.4 Thermal monitoring 

An AD590 temperature sensor is used to monitor the abso- 
lute temperature of the reference load. The reference source 
passively follows the ambient temperature as well as the 
whole instrument. The MPI box and the primary mirror are 
also monitored to estimate possible differential emission in 
the instrument due to temperature gradients. The AD590 is 
a temperature transducer producing an output current pro- 
portional to absolute temperature, suitable for our purposes: 
wide temperature range (-55° C -f- 150° C), high calibration 
accuracy (± 0.5° C) and excellent linearity (± 0.3° C over full 
range) . 
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3.5 Signal modulation techniques 

The peculiarity of CASPER2 is the ability to perform, with 
the same interferometer, 3 different signal modulations: Am- 
plitude Modulation {AM), Fast-Scan (FS) and Phase Mod- 
ulation (PM). The OPD variation along the time for the 
three options is shown in Fig [S] Many fast-scanning FTSs 
are built so that the ZPD position lies close to one edge 
of the moving stage path, in order to maximize the dy- 
namic range of OPD and consequently the spectral resolu- 
tion. Anyway, since one-sided interferograms transform into 
real spectra, no information on the phase is available, al- 
though phase problems do show up as baseline anomalies. 
Two-sided interferograms, on the contrary, transform into 
complex spectra (they have two pieces of information per 
frequency), allowing phase errors to be directly measured as 
a function of frequency. This permits checking for optical 
mis-alignments and other potential instrumental problems 
through the level of asymmetry in the two sides of the in- 
terferograms. Since a high spectral resolution is not neces- 
sary when measuring the continuum level of the atmospheric 
emission, CASPER2 adopts a two-sided interferogram sam- 
pling, with the ZPD located half-way along the moving mir- 
ror path. To reduce the effect of measuring the interfero- 
grams with a limited mechanical path difference, the inter- 
ferogram lobes are weighted with a triangular apodization 
function, at the expense of decreasing the spectral resolution 
to about 8.6 GHz. 

3.5.1 Amplitude Modulation 

The amplitude modulation [AM) technique is fulfilled by 
a chopper wheel modulator placed in front of the cryostat, 
performing a synchronous demodulation of the signal, and a 
step-by-step movement of RM1. The step length is equal to 
100/im ensuring a Nyquist frequency, i/n = 750 GHz, higher 
than 450 GHz, see Table(2] The chopper blades, coated with 
Eccosorb AN72, act as a reference source at ambient temper- 
ature. The two modulated signals feed the two lock-in am- 
plifiers. Since the chopper is kept at room temperature, the 
large brightness temperature gradient between the 2 sources 
implies that the faint signal fluctuations are difficult to de- 
tect. In addition the AM procedure has the disadvantage 
that half of the observation time is spent on the reference 
source. The AM procedure represents a reliable tool useful 
for characterizing the performance of the whole instrument, 
working both with stable sources in laboratory tests and 
instrument calibration over bright sky signals. 

3.5.2 Fast-Scan 

The fast-scan technique consists in sweeping the range of 
available OPDs through a rapid movement of the translat- 
ing stage at constant velocity. The recorded time-domain 
signal is therefore trivially related to the interference pat- 
tern thanks to a simple time/position conversion through 
the stage velocity, and the Fourier analysis yields directly 
the necessary information about the spectrum of the incom- 
ing radiation beam. A carefully selected velocity may shift 
the electrical frequencies of interest away from potentially 
troublesome low-frequency components or line features in 
system noise, allowing for a cleaner reconstruction of the 



optical power in the passband. Moreover, since the scan can 
be repeated an arbitrary number of times, the short integra- 
tion time per unit OPD can be increased to hit the photon 
noise limit with almost no additional effort in the instrument 
setup. On the other hand, apart from mechanical limita- 
tions, an intrinsic upper limit to the value of the translation 
velocity is determined by the time constant of the detectors 
and by the highest frequency in the instrument bandwidth: 
in order to be able to discriminate two consecutive fringes in 
the interference pattern generated by radiation at frequency 
Vmax, they must be scanned in a time interval longer than 
the detector time constant Td, thus determining the limit 
velocity 

Vlim < (15) 

Fast scanning interferometers are usually operated well 
below this limit, and signals are de-convolved from the de- 
tector time-response before further data-processing to avoid 
residual artifacts. This is also the case for CASPER2, where 
the time-constant of the detectors is 7"d=10 ms and the trans- 
lation velocity of the moving roof mirror is set to 0.86 cm 
s _1 (i.e. 1.72 cm -1 lag velocity). Under these conditions, a 
time of 3.8 s is needed to perform one scan. An average of 
several scans on the same source is needed to improve the 
signal to noise ratio of the observation (see Sect. 13.71) . 

A critical source of systematics in fast-scanning FTSs is 
the non-uniformity of the scan velocity: since interferograms 
need to be uniformly sampled in space-domain, rather than 
in time-domain, any fluctuation in velocity at a fixed scan 
rate converts into a local "stretch" of the OPD scale, re- 
sulting in artifacts in the frequency domain, both in line 
and continuum interferometry. In CASPER2, this issue is 
solved by monitoring the stage position through the optical 
sub- interferometer (see Sect. 13. 3j) . whose fringes provide a 
position reference for each data-point. The interferograms 
are highly oversampled (5 kSamples s _1 ), with an electrical 
Nyquist frequency of 2.5 kHz (or an equivalent optical band- 
width of 43.6 THz). After Fourier transformation, the high 
frequency components of instrument noise are discarded and 
the signal in the optically meaningful band is processed to 
extract the sky brightness and pwv information. 

One of the advantages of fast scan interferometry over 
step by step interferometry (or slow scan) is that no chop- 
per is used to modulate the radiation (see Sect. |3"3TT|) . A 
disadvantage in FS interferograms is that slow drifts in the 
intensity of the source can result in variations of the baseline 
of the interferogram which can be of the same frequency as 
modulations from the longest wavelength being measured. 
This affects the performance at low optical frequencies. 

3.5.3 Phase Modulation 

A third signal modulation technique, available with 
CASPER2, is the Phase M o dulati on (PM). This technique, 
proposed by IChamberlainl (|l97ll ). allows the replacement 
of AM and FS implying a modulation of the OPD when 
recording the interferogram. The insensitivity to slow fluc- 
tuations in the intensity of the source plus the full time 
on source observations, make PM very attractive to atmo- 
spheric measurements. 

Operatively modulation is performed by periodically 
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Figure 5. Comparison of variation of OPD along the time for 
FS (green dashed line), AM (red dotted line) and PM (black 
solid line) modulation techniques. 
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Figure 6. The normalized calibration function for CHI (black 
solid line) and for CH2 (red dashed line). 



wobbling RM2 while the detector signals are lock-in de- 
modulated. We can rewrite the Muller matrix for the MPI, 
as in Eq[6] : 



Ah 



t 1 




V 





cos S 


sin S 





sin 8 

-1 

— cos S J 



(16) 



where S — S + 5m with 8m — 2-wM(i) / \, the modula- 
tion of the OPD obtained by periodically oscillating RM2 
with an amplitude A and a frequency vq. In Figure [5] the 
OPD along the time is plotted for all the three observation 
modes. 

The amplitude of the phase modulation has to be care- 
fully chosen to fit the spectral band of the instrument. In our 
case, with a sinusoidal modulation function, A — 0.58c/ 'v max 
that for v m ax = 450 GHz, the highest frequency in the two 
bands, corresponds to 840 fj,m. 

The wobbling frequency is chosen equal to 12 Hz so 
that the effect of low frequency source variations is essen- 
tially eliminated using PM techniques. At the same time, 
by avoiding a dc component the acquisition dynamic range 
is well fitted to the interferogram range. 

RM2 is wobbled by a linear actuator driven by a wave- 
form generator using a feedback loop based on a position 
transducer and a Proportional-Integral-Derivative circuit. A 
linear variable differential transformer (LVDT) (Solartron 
Metrology, model SM/1) is used as position transducer for 
the mirror. The accuracy in oscillation amplitude and fre- 
quency was checked and found to be 2 per cent in amplitude 
and 1 per cent in frequency. The requirements on the stabil- 
ity of the oscillation amplitude were estimated by simulated 
observations. Changes in the amplitude result in a differ- 
ent weighting along the spectrum affecting the inferred pwv, 
as derived by fit of ATM synthetic spectra. The constraint 
on pwv respecting the uncertainty on the RM2 oscillation 
amplitude turns out to be 2 per cent, estimated assuming 
several pwv average values in the range 0.1 mm to 6 mm. 



3.6 Calibration procedures 

Calibrated spectra are derived by employing several consid- 
erations. The two inputs of CASPER2 differ as follows: 

B ini {v) = e atm (v)BB{T atm , v) + e tele (v)BB{T tele , v) (17) 

B ln2 {u) = BB{T ref ,v) (18) 

where e atm and e te ; e are the atmospheric and telescope 
emissivities respectively while BB stands for the specific 
brightness of the atmosphere, the telescope and the reference 
load (the only one having an emissivity equal to 1), each of 
them at the equivalent temperatures T a trn,T te i e and T re f. 
While e a tm is related to the atmospheric opacity, we assume 
tteie equal to 3 • 10 -3 at the frequency of 150 GHz, for an 
aluminum mirror, and changing with the frequency as ^fu 
jBock et all (Il995l )). 

The spectra derived from Eqs. [5] and [5] after baselines 
removal, are related to the spectra of the incoming sources 
as: 

Iou tl (u) = Ri (u)e! (u)An(u) [B ini (u) - B m2 (v)] = 

= Fi (v) [B ini (u) - B in2 (i/)] (19) 

Iout 2 (v) = R 2 (v)e 2 (v)An.(v) [B in » - BimO)] = 

= F 2 {v) [B ln . 2 (i/) - B ini (i/)] (20) 

The calibration function, Fi(y), for the i-channel in- 
cludes the responsivity Ri(y), the spectral efficiency Si(y) 
and the throughput Afl(u). In Fig. [S] the normalized cali- 
bration functions are shown for both channels. 

The throughput AQ(u) is assumed equal for the two in- 
puts, being only dictated by the optical matching between 
the cones and the last cold lens operating as aperture stop. 
A frequency dependence of it can be considered for both 
patterns: starting from a single mode propagation, at the 
cone exit apertures at the lowest frequency, and moving to 
a multi mode approach at 450 GHz. For CASPER2's ap- 
plication this anisotropic response of the two ports is not 
a technical hitch due to the presence of diffuse sources, the 
atmosphere and the sheet of Eccosorb, totally filling both 
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the inp uts. We can easily assume that the telescope effi- 
ciency (|Ruzel (|l966l )) is unitary along the whole spectral 
range counting on an r.m.s. surface error of ~ O.lfim and 
that the wire grids behave like ideal polarizers in our bands. 

The optical matching between the MPI and the two 
sources shows differences due to the distinct optical paths. 
Specifically in the case of in\ , the input port is further trans- 
formed by the primary beam telescope. 

The calibration functions are estimated by filling the ini 
port with a second well modeled source, a cold black body, 
realized with an Eccosorb AN72 sheet thermalized inside a 
liquid nitrogen bath (77 K). Hence the atmospheric spectra 
is deduced by: 

Iat mi {v) = h T} U ? +BB(T ref ,v)~e tele (v)BB(T tele ,is)(2l) 

Possible optical mismatch between the two input ports 
has to be known and taken into account to avoid the conse- 
quent contamination on inferring the atmospheric spectra. 
The correct balance between the two inputs is checked by 
filling even the sky port with a room temperature blackbody 
inserting an Eccosorb AN72 sheet in front of the telescope 
to record a null interferogram. Only in this case a residual 
signal could arise in the acquired spectrum when a far from 
ideal source coupling of the two input ports is present. A po- 
tential temperature gap between the two blackbodies could 
also produce a non-nuJJ interferogram but this however has 
been monitored. 

The knowledge of the null interferogram, or at least 
the upper limit of the ZPD value when the signal-to-noise 
ratio is less then one, allows us to put a constraint on the 
minimum detectable contribution on the pwv content. Long 
acquisition of null interferograms enable us to discriminate 
spectra with a difference of only 0.01 mm of pwv irrespective 
of the pwv content, at least for pwv < 1 mm. 



3.7 The Allan Variance 

The instability of the atmospheric emission in the mm/sub- 
mm spectral region has to be carefully taken into account 
when the observational goal is to achieve frequent and in- 
dependent high signal-to-noise ratio spectra. The time ded- 
icated to performing a single interferogram is mainly dic- 
tated by the detector time response while the timescale 
to average several interferograms is affected by the atmo- 
spheric drifts. In the specific case of the FS technique, to 
avoid loss of observational time, the minimum number of 
spectra that can be averaged has to be constrained achiev- 
ing a high signal-to-noise ratio while sampling a continuous 
changing atmosphere. It is fundamental to determine a char- 
acteristic time to indicate when an instrument is dominated 
only by thermal noise instead of an atmosphere fluctuations 
regime. An appropriate approac h to in f er thi s timescale is 
to estimate the Allan variance (|Allanl l|l966l )). In a wide 
band spectrometer, like CASPER2, it is important to also 
check the timescale similarity for all the frequencies by in- 
vestigating the noise perfo rmance in the measured spectra 
l|Schieder fc Kramerl (|200lT ) V 

Large fluctuations of atmosphere emission are expected 
in correspondence with the 3 "windows" centered at the fre- 
quencies of 150, 270 and 350 GHz; hereafter we refer to these 
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Figure 7. Overplotted atmospheric zenithal spectra recorded at 
MITO on July 16th 2010, starting at 04:03 AM and finishing 
at 04:12 AM. The vertical dotted lines denote the fiducial fre- 
quencies, defined in the text, where the Allan variance has been 
estimated. 



lines with 2, 4 and 5, respectively. On the contrary, the oxy- 
gen band at 118 GHz and the high absorption H2O band at 
183 GHz, quoted as 1 and 3, should appear more stable with 
time. In Fig. [7] for example, 87 spectra of the zenithal atmo- 
spheric brightness recorded by the FS technique are shown 
for the 2 bands of CASPER2. Each spectra is the average 
of a couple of back-and-forth spectra acquired in 6 seconds. 
The vertical dotted lines refer to the examined frequencies. 

The Allan variance is calculated, for the previous 5 ref- 
erence lines, in the following way. The signal, Sf(ti), related 
to the atmospheric emission at the frequency / extracted 
from the FS spectra at time ti, is averaged over variable 
timescales, T, generating the new datasct: 

i+T 

Sf(T,t j ) = - (22) 
i=j+i 

The Allan variance, or the two-sample variance, i s esti- 
mated for the frequency / as in IWiedner et all (|2002l ): 

° a{T) = 2 2- I 2 S f (T, tj ) 

i=2 v 7 

In Fig[8]the Allan variance corresponding to the 5 fre- 
quencies is plotted in the cases of the two CASPER2 bands. 
It is worth noting the expected 1 /T dependence for short av- 
erage times, corresponding to dominant thermal noise, and 
the slope change when atmospheric drifts overcome. This 
behaviour is not satisfied in the case of 1 and 3, where a 
strong and stable emission is present. 

We can estimate, at least for this specific dataset of 
spectra, T ~ 100 s as the more suitable average time. Even 
if the values of the Allan variance increase for all the fre- 
quencies in the high-background Channel 2, affected by a 
larger instrumental noise, the best average time is almost 
the same. We employ this timescale to generate the spectra 
reported in the next Section. 
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Figure 9. Averaged atmospheric zenithal spectra recorded with 
fast scan procedure at MITO on July 16th 2010, 04:03 AM. The 
red line is the best fit obtained with ATM model corresponding to 
pwv = 6.53±0.16 mm for CHI (top panel) and pwv = 6.84±0.16 
mm for CH2 (bottom panel). 



4 PRELIMINARY ATMOSPHERIC SPECTRA 
MEASUREMENTS AT MITO 

In this section we present preliminary measurements of at- 
mospheric spectra performed during the summer campaign 
at MITO from 11th to 19th July 2010. In order to show the 
instrument's capabilities, atmospheric spectra measured by 
the two channels of CASPER2 with fast scan procedure on 
July 16th 2010 are shown in Fig. [9] The uncertainty associ- 
ated with the data points is given by the standard deviation 
of the dataset and it is closely connected to the time on 
which the interferograms average is performed (see Fig. [8}. 
Best fit ATM model corresponds to pwv= 6.53 ± 0.16 mm 
for CHI and pwv= 6.84 ± 0.16 mm for CH2 (red lines in 
Fig. |9]) . The uncertainty on the pwv value was evaluated by 
a random generation of the ATM synthetic spectra within 
the brightness uncertainty range and it corresponds to less 
then 3 per cent. Missing of consistency between data and 
ATM spectra with increasing the frequency, supports the 
necessity to accurately calibrate the dry continuum and the 
H2O pseudocontinu um terms in the simulated atmosphere 
l|Pardo et al.l (120011 )1. The ATM model is currently based 
on T/P typical profiles and the validation of the code is the 
subject of a subsequent paper. 



5 CONCLUSIONS 

A double beam FTS interferometer, installed at the focal 
plane of a 62-cm in diameter telescope, is devoted to mon- 
itoring atmospheric emission spectra in the mm/sub-mm 
band compared to an ambient temperature reference source. 
The MPI is designed to perform three independent signal ac- 
quisition procedures: fast scan, phase modulation and step 
and integrate. Being an ancillary instrument, it is provided 
by an independent altazimuthal mount allowing to point 
the same direction of a telescope dedicated to cosmologi- 
cal millimetre observations from the ground; in the current 
version the 2.6-m in diameter MITO telescope in the Alps. 
The recorded spectra, in the 90-^450 GHz spectral region, 
permi t to validate the res ults of the transfer radiative code, 
ATM (|Pardo et all |2001)), for this site and consequently to 
infer the pwv value as derived by fit with synthetic spec- 
tra or by skydips. The instrument has been characterized 
in laboratory and employed during the observational cam- 
paign at MITO in July 2010. The choice of the best spectra 
integration time is derived by the Allan variance estimate. 
Preliminary spectra recorded by fast scan modulation are 
presented to show the instrument's capabilities. 



ACKNOWLEDGMENTS 

Part of this work was supported by Sapienza Univer- 
sity of Rome with Ateneo 2009 (C26A09FTJ7) and 2011 
(C26A11BYBF). We acknowledge with gratitude F. Mel- 
chiorri's suggestions which were a valid contribution during 
the early phases of this project. Useful proposals were re- 
ceived by G. Pisano, G. Siringo, G. Savini, E.S. Battistelli, 
P. de Bernardis and E. Caca. We also thank the referee for 
the constructive comments which have improved the presen- 
tation of our results and S. Capaldi for the English revision 



of the manuscript. Mechanical and a few optical compo- 
nents of the instrument have been realized at the machine- 
shop of the Department of Physics and at Officina Mecca- 
nica INFN - Sezione di Roma. Preliminary observations were 
performed at MITO with the logistic support provided by 
Sezione IFSI/INAF di Torino, now Osservatorio Astrofisico 
di Torino. 



Atmospheric monitoring with CASPER2 11 



REFERENCES 

Allan D. W., 1966, Proceedings of the IEEE, 54, 221 
Bock J., Parikh M., Fischer M., Lange A., 1995, Appl. Opt., 
34, 4812 

Catalano A. et al., 2004, New A, 10, 79 
Chamberlain J., 1971, Infrared Physics, 11, 25 
Coster A. J., Niell A. E., Solheim F. S., Mendes V. B., 
C. T. P., Buchmann K. P., Upham C. A., 1996, in Pro- 
ceedings of the 9th International Technical Meeting of 
the Satellite Division of The Institute of Navigation (ION 
GPS 1996), Kansas City, MO, pp. 625-634 
De Gregori S., De Petris M., Decina B., Lamagna L., Pardo 
J. R., Petkov B., Tomasi C, Valenziano L., 2012, Monthly 
Notices of the Royal Astronomical Society, 425, 222 
De Petris M. et al., 2005, in EAS Publications Series, 
Vol. 14, EAS Publications Series, Giard M., Casoli F., 
Paletou F., eds., pp. 233-238 
De Petris M. et al., 2007, New A Rev., 51, 368 
Decina B., De Gregori S., De Petris M., Lamagna L., 2010, 
in EAS Publications Series, Vol. 40, EAS Publications Se- 
ries, Spinoglio L., Epchtein N., eds., pp. 107-110 
Gervasi M., de Bernardis P., De Petris M., Masi S., Pisano 

G., 1998, Appl. Opt., 37, 4277 
Matsuo H., Sakamoto A., Matsushita S., 1998, PASJ, 50, 
359 

Matsushita S., Matsuo H., Pardo J. R., Radford S. J. E., 

1999, PASJ, 51, 603 
Paine S., Blundell R., Papa D. C, Barrett J. W., Radford 

S. J. E., 2000, PASP, 112, 108 
Pardo J. R., Cernicharo J., Serabyn E., 2001, IEEE Trans- 
actions on Antennas and Propagation, 49, 1683 
Radford S. J., Holdaway M. A., 1998, in Society of Photo- 
Optical Instrumentation Engineers (SPIE) Conference Se- 
ries, Vol. 3357, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Phillips T. G., ed., 
pp. 486-494 
Ruze J., 1966, IEEE Proceedings, 54, 633 
Schieder R., Kramer C, 2001, A&A, 373, 746 
Serabyn E., Weisstein E. W., 1996, Appl. Opt., 35, 2752 
Wiedner M., Hills R., Pardo-Carrion J., 2002, in Astro- 
nomical Society of the Pacific Conference Series, Vol. 266, 
Astronomical Site Evaluation in the Visible and Radio 
Range, Vernin J., Benkhaldoun Z., Munoz-Tunon C, eds., 
p. 278 

Wiedner M. C, Hills R. E., Carlstrom J. E., Lay O. P., 
2001, ApJ, 553, 1036 



